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A  concept  is  advanced  for  using  the  motion  of  launchers  of  a  free-flight 
launcher/rocket  system  which  is  caused  by  random  imperfections  of  the  rockets 
launched  from  it  to  reduce  the  total  error  caused  by  the  imperfections.  This 
concept  Is  called  ''passive  launcher  control"  because  no  feedback  is  generated 
by  an  active  energy  source  after  an  error  is  sensed;  only  the  feedback 
inherent  in  the  launcher/rocket  interaction  is  used.  Relatively  simple 
launcher  models  with  two  degrees  of  freedom,  pitch  and  yaw,  were  used  in 
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conjunction  with  a  more  detailed,  variable-mass  rocket  model  In  a  digital 
simulation  code  to  obtain  rocket  trajectories  with  and  without  thrust  mis¬ 
alignment  and  dynamic  Imbalance.  Angular  deviations  of  rocket  velocities 
and  linear  deviations  of  the  positions  of  rocket  centers  of  mass  at  burn¬ 
out  were  computed  for  cases  In  which  the  launcher  was  allowed  to  move 
(‘•flexible^  launcher)  and  was  constrained  so  that  it  did  not  rotate  C^rigid” 
launcher)  and  ratios  of  flexible  to  rigid  deviations  were  determined.  Curves 
of  these  "error  ratios"  versus  launcher  frequency  are  presented.  These  show 
that  a  launcher  which  has  a  transverse  moment  of  inertia  about  its  pivot  point 
of  the  same  magnitude  as  that  of  the  centroidal  transverse  moments  of  inertia 
of  the  rockets  launched  from  it  can  be  "tuned"  to  passively  reduce  the 
errors  caused  by  rocket  imperfections.  Additionally,  results  are  presented 
which  show  that  this  reduction  can  be  ah^leved  If  the  imperfections  are 
random.  Hence,  dispersion  of  free-fllght\ rockets  can  be  reduced  by  using 
passive  control  launchers.  Reduction  of  dispersion  by  up  to  eighty  percent 
appears  possible.  \ 
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SECTION  1.  INTRODUCTION 


1.1  Dispersion  of  Er Eliglit  Rockets 

Free-flight,  or  "umont  rol  led, "  loekets  are  very  important  tactical 
weapons  because  they  can  bt'  used  to  rapidly  deliver  large  payloads  with 
good  aci'uracy  and  at  low  I'ost  relative  to  more  sophisticated  weapons.  The 
characteristics  ot  accuracy  and  low  ci>si  arc'  coiillicting  ones.  In  keeping 
costs  low,  concessions  must  he  made  ia  manufacturing  requirements.  This 
results  in  rockets  which  are  more  imperlect  tlian  those  which  could  be  built 
if  cost  was  not  considered.  These  im|)erfect ions  are  random  and  of  many 
types.  They  Include,  for  example', non-st  raiglitness  dynamic  Imbalance, 
static  imbalance,  propellant  imperfections,  thrust  misalignment,  and  geo¬ 
metric  asymmetries.  Since  a  free-flight  rocket  is  not  controlled  after 
it  leaves  the  constraints  of  its  launcher,  random  imperfections  cause  random 
deviations  of  its  actual  flight  patli  from  the  nominal,  or  expected,  flight 
path.  (liven  a  group  of  rockets,  the  summation  of  the  individual  deviations 
is  called  "dispersion." 

Of  course,  there  are  other  tauses  of  dispersion.  For  instance,  one 
immediately  tliinks  of  random  winds.'  However,  thrust  misalignment  and 
dynamic  imbalance  are  two  primary  causes  of  dispersion  of  free-flight  rocket 
In  the  literature  (see,  for  example.  Ref.  2,  pp.  68,  149-156)  the  terms 
"mallaunch"  and  "malaim"  are  used  for  random  transverse  angular  rates  and 
angular  rotations,  respectively,  produced  by  tlie  launcher/rocket  interaction 
Although  changes  in  its  environment  may  cause  a  launcher  to  perform  dif¬ 
ferently  at  different  times,  it  m.ay  be  considered  much  more  deterministic 
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than  the  rockets  launched  from  it.  Hence,  if  the  characteristics  of  a 
particular  launch  are  well  known,  tlie  major  part  of  the  "mallaunch"  and 
"malaim"  associated  with  it  must  he  due  to  imperfect  rockets  which  are 
generally  the  only  energy  source  present. 

If  "mallaunch"  and  "malaim"  are  generated  primarily  by  imperfect 
rockets,  there  is  reason  to  believe  that  -  as  "mal"  denotes  "bad"  -  a 
"goodlaunch"  and  "goodaim"  may  exist,  lie  grounds  for  this  reasoning  are 
that  the  trajectory  of  a  free-fliglit  rocket  is  significantly  affected  by  the 
rocket’s  flight  st.ite  when  it  h'aves  the  influence  of  the  launcher  and  that 
this  state  could,  it  would  appe.-ir,  he  one  which  reduces  the  error  produced 
by  rocket  imperfections  in  free-flight  rather  than  one  which  increases  such 
errors.  The  motivation  for  tliis  investigation  was  provided  by  the  hope  that 
proper  design  of  launchers  will  lead  to  the  replacement  of  "mallaunch"  and 
"malaim"  by  their  antonyms,  whatever  they  may  be. 

1.2  Scope  of  this  Effort 

Previously,  the  connection  between  the  random  part  of  the  motion  of  a 
launcher  and  that  of  a  rocket  launched  from  it  has  been  generally  disregarded. 
Valid  reasons  for  this  are  (1)  that  l.aunchers  are  often  much  more  rmissive 
than  the  rockets  launched  from  them  and  (2)  that  mathematical  simulations  of 
the  motion  of  rockets  and  of  the  motion  of  launchers  have,  for  the  most  part, 
been  pursued  Independently.  This  qualifier  "generally"  is  used  above  because 
there  have  been  Investigations  the  objective  of  which  was  to  design  a  launcher 
which  compensated  for  rocket  Imperfections.  The  PADA^  concept  of  compensat¬ 
ing  for  dynamic  Imbalance  is  a  case  in  point. 


The  goal  of  this  effort  was  to  determine  whether  a  launcher  can  compen¬ 
sate  for  random  imperfections  of  free-flight  rockets.  The  objectives  were 


» 

to  model  the  Inuneher/riicket  system  m.il  Uemal  ical  ly ,  to  use  this  matliemat- 
ical  model  to  determine  qualitatively  how  rocket  imperfections  affect 
launcher  motion  and  to  produce  quantitative  numerical  results  whicli  show  how 
a  launcher  characteristic,  such  as  launcher  Inertia  and  natural  frequency 
of  the  launcher,  affect  its  performenace  ns  a  passive  controller. 

This  goal  was  partially  achieved  during  tl)e  first  phase  of  this  lnves~ 
tigation.’  Tile  results  reported  in  Ret.  4  are  qualitatively  correct,  but 
an  error  in  the  computer  code  which  was  I ound  suhseciuent  to  publication  of 
Ref.  A  renders  the  quantitative  results  tor  launcher/rocket  systems  given 
therein  invalid.  A  later  p;iper  '  pr'csetits  more  accurate  results  which  some¬ 
what  substantiate  tlie  concept  ot  launchers  as  passive  controllers,  but  they 
are  limited  since  only  a  few  cases  were  studied. 

In  the  next  section,  the  concept  of  passive  control  as  applied  to 
rocket  launchers  is  discussed  and  some  practical  considerations  which  are 
warranted  by  physical  laws  ami  intuition  are  stated.  St*ction  3  provides 
descriptions  of  physical  models  of  the  launcher  and  rocket,  as  well  as 
data  used  in  the  numerical  portion  ol  the  investigation. 

The  ratiu’r  extensive  results  contained  in  Section  A  include  plots 
showing  the  effectiveness  of  a  launcher  in  reducing  dispersion  due  to  thrust 
misalignment  and  ilynamlc  imbalance.  I'he  effects  of  launcher  frequency  and 
inertia  are  illustrated  and  the  efiects  of  errors  in  rocket  thrust  profile 
are  quantified  to  some  extent,  tlonclusions  and  recommendations  are  given  in 
Section  5. 

Detailed  descriptions  of  the  mathematical  models  developed  and  used  to 
produce  the  quantitative  results  are  found  In  Appendix  A.  Finally,  details 
of  a  "prototype"  launcher  model  are  given  in  Appendix  B. 


SECTION  2.  PASSIVE  CONTROL 


2.1  The  Concept 

Passive,  as  opposed  to  active ,  control  involves  the  use  of  an  energy 
source,  or  sources,  inherent  within  a  system  and/or  its  environment  and 
"natural"  feedback  to  produce  a  desirable  response  of  the  system  when  a 
disturbance  acts  upon  it.  A  good  ex.ample  of  the  use  of  passive  control  is 
presented  by  gravity-gradient  stabilized  spacecraf  t .  Such  satellites  are 
designed  so  that  tlie  torque  due  to  the  gradient  of  the  earth's  gravitational 
field  causes  the  spacecraft  to  rotate  about  a  centroidal  axis  which  is  per¬ 
pendicular  to  the  satellite's  orbital  plane  once  each  orbit  thereby  pointing 
one  axis  of  the  satellite  toward  the  earth.  In  this  example,  the  "control" 
torque  is  due  to  the  environment  and  the  satellite  configuration. 

Regarding  launclter/ rocket  systems,  the  energy  source  is  the  rocket's 
motor,  and  in  some  cases,  a  spin  motor.  The  feedback  to  the  rocket  is  trans¬ 
ferred  through  the  launcher /rocket  Interface.  All  motion  of  the  system  is 
due  to  the  rocket  motion,  spin  motor,  gravity  and  aerodynamic  forces.  The 
force  due  to  gravity  produces  a  tipping  of  the  launcher  which  is  essentially 
independent  of  the  Imperfections  the  rocket  has  except,  of  course,  the  pro¬ 
pellant  imperfections  which  result  in  different  guidance  times.  Aerodynamic 
forces  Include  those  due  to  the  flow  field  about  rocket  and  launcher  during 
guidance  which  are  very  difficult  to  define,  and  those  on  the  rocket  during 
free-fllght,  which  can  usually  be  determined  quite  well  from  wind  tunnel 
tests. 

The  random  forces  and  moments  of  interest  herein  are  those  due  to  thrust 
misalignment  and  dynamic  Imbalance.  Moreover,  the  thrust  misalignment  is 
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assumed  to  be  angular  misalignment  such  as  that  due  to  a  misaligned  rocket 
nozzle.  Tlie  dynamic  imbalance  is  assumed  to  be  due  to  a  misalignment  of 
the  principal  axis  of  least  inertia  and  the  geometric  symmetry  axis  of  the 
rocket . 

It  can  be  demonstrated* '  tliat  the  flight  errors  due  to  thrust  mis¬ 
alignment  and  dynamic  Imbalance  are  manifested  during  the  time  of  thrusting. 
Tliese  errors  are  of  the  same  type  .is  those  gener.itcd  by  transverse  angular 
rates  of  ttic  rocket  at  tlie  timi*  it  leaves  tlie  launcher.  In  fact,  it  can  be 
shown'  that,  for  a  spinning  rocki-t,  thrust  mis.il  i  gnmcnt  and  dynamic  imbalance 
result  in  transverse  .angular  rates  of  the  rocket  during  the  first  quarter 
of  a  revolution  of  the  rocket  after  It  leaves  the  launcher.  These  angular 
rates  persist  until  damped  out  by  aerodynamic  moments.  Because  of  the 
nature  of  the  errors  due  to  these  two  types  of  rocket  imperfections;  i.e., 
their  equivalence  to  those  due  to  transverse  angular  rates  at  end  of  guidance, 
it  follows  that  if  proper  transverse  angular  rates  of  the  system  are  caused 
by  imperfections  prior  to  the  time  the  rocket  leaves  the  launcher,  then  the 
adverse  effects  of  rocket  imperfections  can  be  reduced.  This  is  the  passive 
control  concept  advanced  herein.  in  simpler  terms,  if  the  imperfections 
cause  the  system  to  rotate  or  to  he  rotating,  in  one  direction  when  the  rocket 
leaves  tlie  launclier  and  they  cause  the  rocket  to  rotate  in  the  opposite 
direction  during  free-flight,  the  It^al  <yt^o_r  will  be  less  than  if  the 
launcher  did  not  move. 

2.2  Practical  Consideration s 

There  are  two  basic  requirements  for  passive  control  by  the  launcher. 
First,  the  imperfections  of  the  rocket  must  cause  the  system  to  rotate. 
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Second,  the  rotation  must  be  in  opposition  to  that  which  the  imperfections 
cause  during  froe-lllght. 

The  first  requirement  dictates  that  either  the  entire  launcher,  or 
the  part  of  the  launcluT  In  contact  witli  the  rocket,  must  be  small 
enough  that  the  imperfections  can  produce  motion  of  it  and  tlie  rocket. 

Hence,  if  the  launcher  (or  the  part  of  it)  is  much  more  massive  than  the 
rocket,  one  would  expect  little  if  any  passive  control  to  be  possible. 

Figure  1  shows  a  concept  of  a  passive  control  launcher  in  which  the  launch 
"tube"  (which  may  not  be  closed,  but  could  be  a  framework)  can  rotate 
transversely  in  response  to  rocket  imperfections.  To  satisfy  the  second 
requirement,  either  the  system  must  rotate  to  point  the  rocket  in  a 
direction  in  opposition  to  that  it  will  rotate  during  free-flight,  or  the 
system  must  acquire  an  angular  rate  which  is  in  opposition  to  that  generated 
subsequent  to  launci).  To  get  sucli  a  system  response,  one  must  design  a 
launcher  (or  part  of  it)  which  Is  supported  so  that  limited  rotation  is 
possible  and,  in  the  case  of  production  of  angular  rates,  has  a  natural 
frequency  which  causes  its  response  to  be  out-of-phase  with  the  post-launch 
response  of  the  rocket. 

An  Important  fact  to  remember  is  that  for  a  spinning  rocket  the  angular 
rate  produced  by  dynamic  Imbalance  or  thrust  misalignment  is  90°  out-of¬ 
phase  with  the  body-fixed  torque  caused  by  the  imperfection.  The  launcher 
response  must  be  determined  accordingly.  Furthermore,  if  the  rocket  is 
not  spinning  when  it  leaves  the  constraints  imposed  by  the  launcher,  the 
duration  of  its  powered  flight  must  not  be  too  much  longer  than  the  time 
it  spent  on  the  launcher:  for,  if  it  Is,  the  angular  rotation  produced  by 
thrust  misalignment  during  free-flight  will  be  much  larger  than  any 


Passive  control  launcher  , 


compensating  rotation  or  angular  rate  that  was  produced  by  the  imperfec¬ 
tion  prior  to  free-f light.  Reference  4  contains  additional  discussion  of 
the  points. 

In  summary,  intuition  and  physical  reasoning  dictate  (1)  a  small 
launcher,  or  a  small  part  of  the  launcher,  which  can  respond  to  the  im¬ 
perfection  and  (2)  proper  phasing  of  the  launcher  and  rocket  motion. 
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SECTION  3.  MODELS 


3.1  General  Comments 

In  this  section  the  models  of  the  launcher  and  rocket  adopted  for 
this  investigation  are  described.  Additional  discussion  of  the  models 
and  the  governing  equations  of  motion  are  provided  in  Appendix  A. 

Although  it  is  somewhat  redundant,  since  similar  comments  are  made 
in  Appendix  A,  it  is  worth  stating  twice  that  the  dynamical  system  of 
vehicle,  launclier  and  rocket(s)  is  a  complicated  one.  Moreover,  the 
variety  of  vehicles,  launcher  configurations  and  types  of  rockets  makes 
the  construction  of  a  model  whicli  would  be  sufficient  for  all  such  systems 
a  very  large  task  at  best.  Although  tlie  development  of  a  fairly  compre¬ 
hensive  model  which  incorporates  major  characteristics  of  launcher/rocket 
systems  seems  justified,  it  is  generally  more  efficient  for  certain 
individuals  to  develop  the  skills  needed  to  construct  mathematical  models 
of  launcher/rocket  systems  and  to  apply  these  to  the  particular  concept 
under  study. 

Often,  simple  models  of  a  dynamical  system  yield  extremely  useful 
results.  What  must  be  modeled  are  the  salient  characteristics  of  the  system, 
not  all  the  details.  The  models  described  in  what  follows  are  of  this 
type.  Tiiey  are  relatively  simple  (especially  t)ie  launcher  model,  and  yet 
provide,  we  believe,  a  good  idea  of  the  dynamical  behavior  of  the  class 
of  systems  considered. 

3.2  Launcher  Model 

The  "launcher"  model  is  actually  a  model  of  either  a  small  light-weight 
single-round  launcher,  or  what  we  term  a  "sublauncher."  Due  to  the  require- 
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meiii  tliat  tliero  must  be  a  response  to  the  relatively  small  forces  and/or 
either  the  launclier  must  be  small  or  there  must  be  a  relatively  small 
part  which  responds.  It  is  assumed  that  In  the  latter  case  the  motion 
of  the  small  part,  or  sublauncher,  due  to  rocket  imperfections  does  not 
cause  significant  motion  of  the  main  portion  of  the  launcher.  Regardless 
of  whether  it  is  a  launcher  or  a  sublauncher,  we  shall  refer  to  it  as 
"the  launcher." 

The  model  is  a  tube,  or  framework,  which  has  two  degrees  of  freedom, 
pitch  and  yaw.  Two  examples  are  shown  in  Fig.  2.  Both  of  these  are  con¬ 
figured  so  that  they  rotate  about  a  point  on  their  longitudinal  axes. 

The  rotation  is  constrained  by  "leaf  springs."  In  both,  helical  rails  are 
assumed  to  provide  a  specific  angular  rotation  of  the  rocket  about  its 
longitudinal,  or  roll,  axis  during  guidance.  It  ahoul-d  be  emphaai'^.ed 
that  helioaL  vails  are  not  required.  However,  specific  roll  angles  were 
considered  in  obtaining  results  presented  infra  and  helical  rails  are  mt!ans 
of  providing  these.  The  results  are  valid  if  other  methods  for  producing 
spin,  such  as  spin  vanes,  are  used  as  long  as  the  angle  of  rotation  can  be 
controlled  closely. 

The  launcher  is  assumed  to  have  axisyiranetric  inertia  characteristics 
and  to  have  equal  natural  frequencies  in  pitch  and  yaw.  Viscous  damping  is 
assumed . 

The  launcher  is  completely  specified  by  M,  its  mass;  1^^,  its  trans¬ 
verse  moment  of  inertia  about  its  center  of  mass  C^;|r|,  the  distance  from 
C  to  the  pivot  point  0  (see  Appendix  A,  Fig.  A2);  w  ,  its  undamped  natural 

La  wil 

frequency;  and  damping  ratio. 
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Tlie  launclier's  moment  of  Inertia  about  its  longitudinal  axis  is  im¬ 
material  because  it  does  not  rotate  about  that  axis.  A  value  of  R<0 
implies  that  the  pivot  point  is  forward  of  the  launcher's  center  of  mass. 

[n  the  following  section,  several  different  pivot  locations  are  con¬ 
sidered.  The  launcher  characteristics  are  given  in  Table  1.  To  simplify 
matters,  in  all  cases  except  those  for  the  prototype  launcher,  R=0  and  the 

pivot  location  is  determined  by  the  vector  s  ,  from  the  launcher  ccntei'  of 

-a 

mass  to  point  a  on  the  rocket. 

3.3  Rocket  Model 

The  rocket  model  is  necessarily  more  detailed  than  that  for  the  launcher. 
The  main  reasons  for  this  are  that  (1)  the  mass  of  the  rocket  varies,  (2)  the 
rocket  is  subjected  to  aerodynamic  forces  which  must  be  described  in  some 
detail,  and  (3)  the  rocket  is  the  part  of  the  system  which  has  imperfections. 

The  rocket  is  modeled  as  a  variable-mass  body  which  is  rigid  except  for 
the  fluid  wltliin  which  is  generated  by  burning  solid  propellant.  Aside 
from  the  variation  of  the  rocket's  mass  and  moment  of  inertia  with  time, 
the  variability  of  its  mass  is  assumed  to  cause  only  thrust  and  a  jet  damp¬ 
ing  torque  during  free-f light.  Other ^  reaction  forces  and  torques  due  to 
internal  combustion  and  fluid  flow  are  assumed  negligible. 

The  propellant  charge  is  modeled  as  a  right-circular  cylinder  with  a 
cylindrical  port  for  the  purposes  of  determining  the  moments  of  inertia  of 
the  propellant  at  time  t  and  obtaining  an  expression  for  the  jet  damping 
torque.  The  Internal  radius  changes  so  that  the  mass  of  propellant  at 
time  t  agrees  with  that  computed  from  the  thrust  equation  (pressure  term 

omitted)  dm/dt  =  -F„,(t)/V  ,  where  m  is  the  rocket  mass;  F  is  the  thrust 
1  L  I 
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magnitude,  whicli  is  tabulated;  and  is  the  characteristic  exhaust  velocitv 
(speed)  of  the  propulsive  subsystem.  The  thrust  profile,  F.j,(t),  used  in  this 
investigation  as  well  as  rocket  mass,  moments  of  inertia  and  dimensions  are 
given  in  Table  2.  They  are  typical  data  compiled  from  data  for  an  experi¬ 
mental  rocket.'^ 

The  aerodynamic  forces  and  moments  are  modeled  fairly  simply.  As  ex¬ 
plained  in  Appendix  A,  they  depend  on  the  flight  Mach  numbers,  aerodynamic 
angles,  the  angular  velocity  of  the  rocket  and  its  speed.  No  wind  effects 
are  considered.  The  aerodynamic  coefficients  used  are  listed  in  Table  3. 


I  I 


'lab  le 

3  Roc  ke  t 

aerodynamic  characteri.stics 

Re  ference 

Area  (m' ] 

I;  0.0192 

3  Rcierciice  l.cnglh  (m)  :  0.157 

Mach  No. 

c; 

X  (m  from  nose) 

N 

m 

cp 

'  1 

q 

0.0 

4.41 

-1460 

1.876 

0.4 

4.24 

-1540 

1.722 

0,6 

4.12 

-1610 

1.628 

0.8 

3.95 

- 1 7  30 

1.457 

0.9 

3.80 

-1780 

1.408 

1  .0 

3.58 

-1780 

1  .  384 

1.1 

3.84 

-17  30 

1  .378 

1.2 

4.07 

-1660 

1.378 

1.4 

4.41 

-1560 

1.423 

1.6 

4.61 

-1460 

1.469 

1.0 

4.70 

-1380 

1.512 

2.0 

4.76 

-1300 

1.533 

2.2 

4.71 

-1230 

1.548 

2.4 

4.67 

-1170 

1.533 

2.6 

4.54 

-1120 

1.487 

3.2 

4.41 

-1010 

1.298 

Mach  No. 

0,00  0.425 

0.80  0.328 

0.90  0.310 

0.95  0,305 

1.00  0.330 

1.10  0.403 

1.20  0.381 

1.28  0.373 

1.50  0.372 

2.00  0.340 

2.50  0.299 

3.00  0.262 

3.50  0.231 

4.00  0.205 


SKCTION  4.  PASSIVE  CONTROI,  RKSIM.TS 


4.1  Definition  of  Imperfection 

The  effects  of  two  types  of  Imperfection,  dynamic  Imbalance  and  mis¬ 
alignment  of  the  thrust  force,  are  considered  in  this  section.  The  dynamic 
imbalance  is  defined  by  li.,  and  (see  Appendix  A),  whicli  are  small  angles 
generated  by  rotating  from  a  centroidal  principal  system,  Cx*y*z*,  to  a 
rocket-fixed  system,  Cxyz,  wliich  has  its  x-axis  parallel  to  the  symmetry 
axis  of  the  launcher.  Figure  3  illustrates  the  angle  linear  mis¬ 

alignment  of  the  thrust  force  with  respect  to  the  symmetry  axis  of  the  rocket 
is  defined  by  the  vector  from  point  "a"  on  the  symmetry  axis  of  the 

rocket  to  a  point  an  the  line  of  action  of  the  thrust.  Angular  thrust  mis¬ 
alignment  is  specified  by  the  small  angles  i  and  i  which  are,  effectivelv, 

y  z 

rotations  about  the  y  and  z  rocket-fixed  axes,  respectively.  In  this  report, 

■.  has  only  an  x-component .  Hence,  all  the  misalignment  is  due  to  >  and  a  . 
-i  y  z 

Figure  4  illustrates  the  angle 

For  all  the  dynamic  imbalance  results  presented  lierein,  except  those  in 

subsection  4.9,  =  0.001  and  =  0.  Also,  with  tlie  above  exception,  for 

all  the  thrust  misalignment  cases,  t  =  -0.001  and  i  =  0. 

y  z 

4.2  Definition  of  Flight  Error.s 

Since  dynamic  imbalance  and  thrust  misalignment  cause  changes  in  the 
direction  of  flight  of  free  rockets  during  the  powered  portion  of  free- 
f light,  they  result  in  two  types  of  errors  at  burnout.  One  is  in  the 
direction  of  the  velocity  of  an  affected  rocket.  The  second  is  in  the  dis¬ 
placement  of  its  center  of  mass  from  the  desired  trajectory.  The  error 
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in  flight  direction  is  the  most  significant  if  the  unpowered  phase  is  much 

longer  than  the  powered  phase,  since  it  propagates  as  a  displacement  erriir. 

The  direction  of  flight  at  any  instant  is  defined  by  the  angles  >  and 

,  whicti  orient  the  rocket's  velocity  vector  (see  Fig.  5).  The  angle  i 
w 

is  measured  from  the  liorizontal  plane  in  the  vertical  plane  which  contains 

the  rocket's  center  of  mass  and  the  velocity  V^.  The  angle  j  is  used  to 

'  -C  w 

define  the  orientations  of  tlie  vertical  plane  with  respect  to  some  reference 
line,  perhaps  the  horizontal  projection  of  the  desired  line  of  flight. 

The  position  of  the  center  of  mass  of  a  rocket  is  given  by  the  coordi¬ 
nates  X,  Y  and  7.,  where  X  is  nominally  measured  along  the  desired  direction 
of  flight,  Z  is  measured  positive  down  (a  flat  earth  is  assumed)  and  Y  is 
measured  positive  to  the  right  of  the  desired  flight  direction. 

In  order  to  define  errors  using  the  angles  )  and  ij.i  and  the  displace- 

w 

ments  Y  and  Z,  we  must  have  a  standard.  The  obvious  one  is  the  unperturbed 
or  nominal  trajectory;  i.e.,  the  trajectory  flown  by  a  perfect  rocket.  By 
letting  a  subscript  N  denote  a  nominal  value  and  P  a  perturbed  value  of  a 
state  variable,  we  have  the  following  errors  as  functions  of  time: 

A>(t)  =  Yp(t)  - 

Ai^  =  <p  (t)  -  Ip  (t) 

W  V  w 

P  N 

ay  =  Yp(t)  -  Yj^(t) 

AZ  =  Zp(t)  -  Zj^(t) 

The  total  angular  deviation  or  error,  is  /Ky^+Kip^  and  the  total  transverse 
linear  error  is  /AY^+AZ^  .  The  purpose  of  passive  control  is  to  reduce 
these  two  errors  at  burnout  regardless  of  the  random  nature  of  their  cause 
and  hence  reduce  the  total  dispersion  of  the  rockets  launched. 
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4.3  Determination  of  Launcher  Performance 


Some  method  must  be  established  for  determining  how  well  a  given 
launcher  performs  in  its  passive  control  role.  The  following  method  has 
been  adopted; 

1.  A  specific  launcher  configuration  is  chosen. 

2.  A  specific  type  of  Imperfection  is  chosen. 

3.  A  nominal  trajectory  is  generated  by  numerically  integrating 

the  equations  of  motion  of  the  system  with  no  imperfections 

and  a  rigid  launcher;  i.e.,  w  =  until  burnout  (t=t,  ). 

■  Ln  bo 

4.  A  perturbed  trajectory  is  generated  by  numerically  integrating 
the  equations  of  motion  of  the  system  with  the  specific  imper¬ 
fection  and  a  rigid  launcher  until  burnout. 

5.  Step  3  is  repeated  with  lOj  equal  to  a  finite  number;  i.e.,  a 
flexible  launcher. 

6.  Step  4  is  repeated  with  equal  to  the  value  in  step  5. 

7.  The  difference  in  the  flight  path  angles  Y,  and  ijJ  (Ay  and 

bo  w, 

bo 

AiJj  ,  respectively)  and  of  the  lateral  displacements  Y,  and 
w  bo 

bo 

Z,  (Ay,  and  AZ.  ,  respectively)  are  calculated  from  the 
bo  bo  bo 

results  of  steps  3  and  4  (rigid  launcher  errors)  and  from  the 
results  of  steps  5  and  6  (flexible  launcher  errors). 


8.  The  ratios  of  /Ay^.  +  and  /AY^  +  AZ/;  for  the  flexible 

bo  w,  bo  bo 

bo 

launcher  to  the  same  quantities  for  the  rigid  launcher  are 
computed . 

9.  Steps  6  through  8  are  repeated  with  other  launcher  frequencies 
to  generate  launcher  response  curves. 

10.  The  entire  process,  steps  1  through  9,  is  repeated  for  a  different 
launcher  configuration  and/or  imperfection. 

The  foregoing  procedure  in  no  way  guarantees  that  an  optimal  launcher 

configuration  will  be  found.  However,  if  enough  of  the  parameter  space 
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spannod  liy 
angle  (]>  is 
cliallenge , 


Llie  pliys  i  la  1  I  y  p<issil>li-  values  of  1^  ,  ,  K  and  rocket  rotation 

searched  at  least  we  may  *)btain  results  wliiclt  will  support,  or 
the  passive  control  concept. 


4.4  Passive  Control  of  Dispersion  due  to  Thrust  Misalignment 

Obviously,  only  a  finite  number  I'f  parameter  sets  could  be  considered 
in  the  finite  amount  of  time  available  for  this  investigation.  This  is 
doubly  true  because  of  tlie  amount  of  computer  time  required  to  produce 
nominal  and  perturbed  trajectories.  in  ail,  six  sets  of  thrust  misalignment 
error  ratio  versus  launcher  natural  frequency  curves  were  produced  for  a 
f ixed  launcher  transverse  moment  of  inertia,  =  54.24  kg-m‘ .  This  moment 
of  inertia  was  fixed,  so  that  the  effects  of  geometry;  i.e.,  launcher  support 
location  and  rocket  rotation  angle  could  be  studied. 

The  first  set  of  curves  for  an  aft  launcher  support  location  and  a  non¬ 
tip-off  guidance  length  of  1.524  m  are  sitown  in  Fig.  6.  Three  different 
rotation  angles  were  used  in  producing  the  curves.  The  first  thing  to  note 

about  this  set  is  that  for  an  aft  support  location  (s  =0;  see  Appendix  A) 

-a 

even  the  "best"  launclier  reduces  the  error  by  less  than  10  percent.  Since 

we  have  plotted  error  ratios  for  equal  rotation  angles  at  end-of-gu idance 

(liOG)  tlie  curves  may  he  compared  even  Lhougli  lor  the  larger  the 

absolute  error  would  be  less  because  tlie  rocket  is  spinning  more  rapidly  at 

KOG.  Note  also  that  the  curves  for  +Ay,  and  +AZ,^  error  ratios 

w.  bo  bo  bo 

b 

are  essentially  the  same. 

When  the  support  is  moved  to  the  front,  the  launcher  becomes  more  sen¬ 
sitive  to  tlie  imperfection.  Kvidence  of  this  sensitivity  can  be  seen  in 
Fig.  7.  Note  that  at  a  launcher  frequency  of  around  18  Hz,  the  error  ratio 
for  ((>„„„  =  ff/2  Is  only  about  0.2.  Therefore,  an  error  reduction  of  about 

EOvf 

80  percent  is  predicted  for  such  a  launcher. 


0  5  10  1$  20  25  30 
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Fig.  7  Error  ratios  -  thrust  misalignment,  forward  pivot, 
non-t Ip-of  f . 
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Tip-off  launcher  results  are  given  in  Fig.  8.  Curves  for  =  tt/2 

and  are  given.  The  tip-off  distance  is  0.3048  ra  and  EOG  occurs  at  the 

beginning  of  tip-off.  The  launcher  pivot  point  is  again  at  the  rear. 

Because  the  reduction  for  the  case  is  greater  than  10  percent  for 

EO(j 

all  launcher  frequencies,  the  tip-off  apparently  is  beneficial. 

When  the  tip-off  distance  is  increased  (see  Fig.  9)  to  1.524  m,  rocket 

rotation  so  that  =  i'/2  is  more  beneficial  and  the  error  ratio  is 

EOG 

drastically  reduced  to  less  than  0.5  at  oj  12  Hz. 

l.n 

The  front-pivot,  tip-off  results  shown  in  Figs.  10  and  11  are  of 
interest  also.  The  0.3048  m  tip-off  curves  look  much  like  the  corresponding 
non-tip-off  curves  of  Fig.  6.  However,  tl'.e  1.524  m  tip-off  curves  are  very 
different  from  their  non-tip-off  counterparts  and  indicate  that  tip-off  is 
detrimental  if  the  launcher  is  pivoted  near  its  front  end. 

In  addition  to  the  constant-moment-of- inertia  results  given  in  Figs.  6 
throgh  11,  results  were  obtained  for  the  "prototype''  launcher  which  is  de¬ 
scribed  in  Appendix  B.  This  launcher  Is  basically  a  steel  tube.  Its 
centroldal  transverse  moment  of  Inertia  is  61.17  kg-m^ .  Only  non-tip-off 
results  were  obtained  using  this  model.  These  are  shown  in  Figs.  12,  13 
and  14  which  correspond  to  aft  ([x  j  =-3.48  m) ,  mid  (  jx^  |  =-1.454  m)  and 

forward  ( I  “  0.07  m)  point  locations,  respectively.  The  launcher  pivoted 
a 

at  its  front  is  obviously  superior.  Especially  encouraging  is  that  if  f  =2  , 

E0(' 

a  launcher  with  =  7  Hz  is  a  very  good  passive  controller.  The  launcher 
supported  at  its  center  of  mass  also  acts  as  a  passive  controller,  but  is 
less  effective.  Furthermore,  the  launcher  which  is  pivoted  near  its  aft 
end  does  not  respond  as  well  as  the  other  two.  A  maximum  error  reduction 
of  only  about  12  percent  was  achieved  with  the  aft  pivot.  This  is  at  least 
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Fig.  ii  Error  ratios  -  thrust  misalignment,  torward  pivot,  large 
tip-off  distance. 
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Fig.  13  Prototype  launcher  error  ratios  -  thrust  misalignment, 
mld-polnt  pivot. 
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partially  due  to  the  increase  of  the  inertia  of  the  system  during  launch 
if  the  pivot  point  is  at  the  rear  of  the  launcher. 

4.5  Passive  Control  of  Dispersion  due  to  Dynamic  Imbalance 

An  axisynunetr ic  rigid  body  wliich  is  constrained  to  spin  about  an  axis 
which  is  not  aligned  with  its  axis  of  minimum  or  maximum  centroidal  moment 
of  inertia  and  then  released  and  allowed  to  rotate  freely  nutates  and  pre- 
cesses.”  The  nutation  has  a  period  close  to  the  spin  period,  but  if  the 
body  is  a  very  slender  body  of  revolution,  the  precessional  period  is  much 
longer  than  the  spin  period.  In  a  rocket  which  is  dynamically  unbalanced 
the  long  period  precession,  coupled  with  the  thrust  force,  causes  the 
rocket  to  depart  from  its  nominal  trajectory. 

Dynamic  imbalance  of  a  rocket  which  is  rotating  on  its  launcher  pro¬ 
duces  a  torque  on  the  launcher.  If  the  launclier  is  sufficiently  responsive 
to  this  torque,  beneficial  motion  of  the  system  may  be  achieved. 

Ten  sets  of  launcher  response  curves  for  dynamic  imbalance  and  fixed 
launcher  moment  of  Inertia  are  presented  in  this  subsection.  Tliey  are 
divided  into  non-tip-off  and  tip-off  groups. 

Non-Tip-Off  Results 

Figure  15  gives  the  response  corves  for  ^  launcher 

E()t» 

with  an  aft  pivot  point.  A  great  deal  of  passive  control  is  not  evident. 
However,  the  ratios  are  always  less  than  unity.  Thus,  the  launcher's 
flexibility  is  beneficial  regardless  of  its  frequency. 

A  mid-point  pivot  was  also  considered.  From  Fig.  16  it  appears  that 
the  mid-point  pivot  is  considerably  better  than  an  aft  pivot.  However, 

r 

the  angles  of  rotation  used  to  generate  Fig.  16  are 


=  371/ 2  and  271. 
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For  direct  comparison  we  need  Fig.  17,  tlie  dashed  curve  of  which  can  be 


compared  with  the  solid  curve  of  Fig.  16.  Upon  making  this  comparison,  it 
is  clear  that  a  mid-point  pivot  is  better  than  an  aft  pivot. 

Use  of  a  forward  pivot  also  leads  to  positive  results,  as  shown  by 
Fig.  18.  At  low  launcher  frequencies,  the  error  ratio  is  much  less  than 
unity.  The  minimum  is  about  0.3. 

Tip-Off  Results 

The  results  for  an  aft  pivot  and  a  tip-off  distance  of  0.3048  m  (see 
Fig.  19)  are  about  as  good  as  the  corresponding  non-tip-off  results.  How¬ 
ever,  the  1.524  m  tip-off  distance  curves  shown  in  Fig.  20  never  pass  below 
the  break-even  value  of  1.  The  same  statements  hold  regarding  the  mid¬ 
point  pivot  tip-off  results  shown  in  Figs.  21  and  22. 

The  forward-pivot  results  for  0.3048  m  tip-off  presented  in  Fig.  23 
are  not  nearly  as  good  as  those  shown  in  Fig.  18  (non-tip-off)  except  at 

very  high  frequencies  and  =  n/2.  The  1.524  m  tip-off  results  (see 

EOC» 

Fig.  24)  are  better.  For  =  ''/2  and  ui,  =  2  Hz,  the  error  ratio  is 

EOC  Ln 

about  0.3. 

The  "prototype"  launcher  results  presented  in  Figs.  25,  26  and  27 
clearly  reveal  the  sensitivity  of  a  launcher  with  a  mid-  or  forward-pivot 
point  to  dynamic  imbalance.  These  curves  arci  all  for  non-tip-off  launchers. 

4.6  Passive  Control  of  Dispersion  due  to  Combinations  of  Imperfections 
Generally,  both  thrust  misalignment  and  dynamic  imbalance  will  be 
present.  Thus,  to  be  successful,  the  launcher  must  passively  control  random 
combinations  of  these  two  types  of  Imperfections. 

The  "prototype"  launcher  model  was  used  to  determine  launcher  per¬ 
formance  in  the  presence  of  thrust  misalignment  ('i  =  -0.001,  u  =  0)  and 

y  z 


I  «  lilt 


0 

5  10 

IS  20 

Of  (HZI 

25 

30 

Ln 

Fig,  18 

Error  ratios 
tip-off. 

-  dynamic  imbalance, 

forward 

pivot,  non- 

I 


47 


dynamic  imbalance  ~  0.001  and  ii.^  =  0).  The  curves  in  Figs.  28,  29  and 
30  indicate  that  positive  control  is  possible.  A  mid-point  pivot  is  best 
for  a  range  of  values  of 

LOG 

Since  dynamic  imbalance  effects  Increase  with  spin  rate,  while  thrust 
misalignment  effects  decrease,  the  results  for  =  it  and  2:  were 

principally  produced  by  dynamic  imbalance.  It  is  not  surprising  there¬ 
fore  that  ttiey  are  very  similar  to  those  for  dynamic  imbalance  only. 

4.7  Effect  of  Launcher  Inertia 

Some  indication  of  the  effect  of  launcher  moment  of  inertia  can  be  de¬ 
rived  from  the  "prototype"  launcher  model  results.  However,  to  delineate 
a  range  of  moments  of  inertia  for  which  passive  control  may  be  feasible, 
the  best  cases  from  several  of  the  "constant-inertia"  results  error 
ratios  were  chosen  and  error  ratios  were  determined  for  various  values  of 
I, /!„»  where  I  is  the  initial  value  of  the  rocket’s  transverse  moment  of 
inertia.  As  can  be  seen  from  Fig.  31,  in  which  results  for  several  thrust 
misalignment  cases  are  presented,  the  system  performance  vapidly  deterior¬ 
ates  as  the  ratio  I, /!„  is  increased.  Figure  31  is  for  non-tip-off  launch- 
L  K 

ers.  For  a  tip-off  launch  distance  of  0.3048  m,  the  ratio  I,/I„  can  be 

L  R 

larger.  Note,  however,  the  high  launcher  frequencies  which  were  used  to 
obtain  Fig.  32. 

Figure  33,  for  non-tip-off  launchers  and  dynamic  Imbalance,  provides 
more  encouragement  as  far  as  launcher  size  is  concerned.  If  a  mid-point 
pivot  is  used,  it  appears  that  a  relatively  large  launcher  will  still 
respond  significantly  to  dynamic  Imbalance. 
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Fig.  30  Prototype  launcher  error  ratios  -  combination  oi 
Imperfections,  forward  pivot,  non-tip-off. 
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fig.  32  fc’flecL  ut  launcher  inertia  -  thrust  misalignment, 
forward  support,  small  tip-off  distance. 
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4.8  Effect  of  Thrust  Profile  Variations 

A  brief  look  was  taken  at  the  effects  of  thrust  profile  variations  on 
launcher  performance.  A  simple  tlirust  profile  (see  Fip.  34)  was  assumed 
and  the  thrust  rise  time  *^2”*‘l*  varied  while  keeping  the  total  impulse 

constant.  The  differences  in  error  ratios  from  a  nominal  profile  were 
determined  for  various  launcher  frequencies.  Figure  35  contains  plots  of 
the  percent  error  in  the  error  ratios  based  on  =  0.06  sec  as  the  nominal 
value.  For  all  cases,  t^^  =  0.005  sec  was  used.  Dynamic  imbalance  was  the 
imperfection.  The  launcher  was  the  basic  one  and  had  a  forward  pivot  point. 
The  launcher  frequency  was  varied  as  indicated. 

Note  that  even  relatively  large  changes  in  tiie  thrust  rise  time  result 
in  only  small  percentage  changes  in  the  error  ratios.  Tims,  the  passive 
control  capability  of  the  launcher  should  not  be  significantly  affected  by 
small  random  clianges  in  thrust  profile  from  rocket  to  rocket. 

4.9  F.ffects  of  Randomness  of  Imperfections 

Since  the  results  presented  in  the  previous  subsections  wore  obtained 

by  using  one  or  more  specific  imperfections,  there  is  obviously  a  question 

as  to  whether  they  are  valid  for  random  imperfections.  One  way  to  check  to 

see  if  they  are  is  to  vary  the  rocket-fixed  orientation  of  the  imperfectionCs) 

through  2ti  and  produce  error  "ellipses."  These  should,  if  the  magnitude  of 

the  Imperfection  is  not  varied  and  there  are  no  biases,  be  circles  centered 

at  the  origin  of  Ay,  vs.  Aijj  and  Ay.  vs.  AY 

bo  w,  bo  bo 

bo 

Figure  36a  shows  typical  error  curves.  The  solid  curves  w-re  obtained 
for  a  thrust  misalignment  case  for  which: 

^  =  -0.001,  =  0,  =  2,1,  =  8  llz. 


the  pivot  point  was  at  tlio  I  rotU  I'l  tin-  laiiiu  lu'i  .  tlu'ic  was  lu-  i  Ip  .Ml  an. I 
I, /l„  =  0.5.  The  dashed  curves  are  for  a  rigid  launcher  and  a  rocket 

r-  - 

with  the  same  thrust  misalignment  and  value  of  (Ji  „•  A  large  reduction 

EOG 

in  total  error  is  clearly  evident. 

The  plots  in  Fig.  Thb  ;ire  for  the  same  case  except  I, /!„  “  6.  The 

L  R 

"overshoot"  exhibited  in  Flp.  51  appears  here  in  the  form  of  flexible 
launcher  error  "ellipses"  which  are  larger  than  their  rigid  launcher 
counterparts . 

Error  curves  for  a  dynamic  imbalance  case  are  presented  in  Fig.  36. 

For  all  the  curves,  =  0.001,  ii  =  0,  4'  ~  and  a  tip-off  distance 

2.  J  EOG 

of  0.3048  m.  The  dashed  curves  are  again  for  a  rigid  launcher.  To  obtain 

the  solid  curves,  a  launcher  frequency  of  6  Hz  and  a  forward  pivot  point 

were  used.  The  dashed  curves  in  Fig.  37a  and  Fig.  37b  are  obtained  using 

I  /I  =  0.5  and  I, /!„  =  8,  respectively  The  large  reduction  in  dispersion 

L  K  L  K 

for  small  I, /l„  and  the  smaller  reduction  for  l./I,  are  consistent  with 
L  R  b  R 

previously  presented  results. 

There  is  no  observable  bias  in  the  curves  presented  in  Figs.  36  and 
37  and  they  clearly  show  that  the  body-fixed  orientation  of  the  imperfec¬ 
tion  does  not  affect  the  launcher  performance  as  a  passive  controller. 

The  effect  of  varying  the  magnitude  of  the  imperfection  is  simply  a 
change  in  the  sizes  of  both  the  rigid-launcher  and  flexible-launcher  error 
curves.  Hence,  the  randomness  of  the  imperfections  does  not  affect  the 
launcher’s  performance  as  a  passive  controller. 


SECTION  5.  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Conclusions 

The  goal  of  this  investigation  was  to  determine  whether  a  launcher 
can  compensate  for  random  imperfections  in  free-fllght  rockets,  thereby 
reducing  dispersion.  This  goal  was  reached  to  the  extent  that  reduction 
of  flight  errors  for  particular  launcher  configurations  was  demonstrated 
by  numerically  solving  the  equations  of  motion  of  the  launcher /rocket 
system  model.  The  reduction  in  error  was  found  to  be  as  large  as  eighty 
percent  for  launcher  natural  frequencies  and  moments  of  inertia  which 
appear  to  be  physically  realizable. 

The  objectives  of  modeling  the  launcher /rocket  system  and  producing 
quantitative  results  were  attained.  The  quantitative  results  indicate 
that: 

1.  A  forward  pivot  point  for  the  launcher  in  general  results  in 
more  significant  response,  especially  for  thrust  misalignment; 
while  either  forward,  or  mid-point,  pivot  is  best  if  dynamic 
imbalance  is  present. 

2.  The  natural  frequencies  of  launchers  which  produce  significant 
reduction  in  dispersion  when  thrust  misalignment  is  present  are, 
for  some  launcher  configurations,  low  enough  to  be  physically 
realized;  while  if  dynamic  imbalance  is  present  it  appears  that 
any  launcher  flexibility  is  beneficial. 

3.  Tip-off  launchers  produce  larger  reductions  of  error  when  the 
imperfection  is  thrust  misalignment,  but  do  not  do  so  if 
dynamic  imbalance  is  the  imperfection. 
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4.  The  effect  of  increasing  launcher  inertia  is  (as  intuition 
implies)  to  reduce  the  capability  of  the  launcher  to  respond 
to  the  force  and/or  torque  caused  by  the  imperfection. 

5.  The  effects  of  small  variations  in  thrust  rise  time  are  not 
significant  as  far  as  launcher  performance  is  concerned. 

6.  The  randomness  of  imperfections  does  not  affect  launcher  per¬ 
formance  . 

5.2  Recommendations 

On  the  basis  of  the  investigation,  we  make  the  following  recommenda¬ 
tions  : 

1.  The  results  obtained  herein  should  be  verified  by  using  analog 
simulations  and/or  other  digital  computer  codes. 

2.  A  simple  experimental  launcher  should  be  constructed  and  tests 
made  under  carefully  controlled  conditions  to  experimentally 
substantiate  the  passive  control  concept. 

3.  Additional  analytical  work  should  be  undertaken  to  determine 
the  theoretical  performance  of  a  passive  control  launcher  which 
has  more  degrees  of  freedom  and  contains  a  "sublauncher"  (see 
Fig.  1).  The  flexibility  of  the  rocket  might  also  be  modeled  to 
ascertain  its  effect  on  launcher  performance. 
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APPENDIX  A 


MATHEMATICAL  MODELS 


A.l  Introduction 

This  appendix  contains  a  description  of  the  mathematical  models  used 
to  obtain  the  results  described  in  the  body  of  this  report.  Vector /matrix 
notation  used  in  the  derivation  and  manipulation  of  equations  of  motion 
is  first  defined.  Then, the  physical  models  for  the  rocket  and  the  launcher 
are  described.  Finally,  equations  of  motion  (mathematical  model)  are 
derived  and  cast  into  the  matrix  forms  which  have  been  programmed  for 
numerical  solution. 

A. 2  Vector /Matrix  Notation 

The  notation  used  for  vector  quantities  such  as  position,  velocity 
and  angular  velocity  is  typified  by 

/V  ^  /S 

-  "  ‘*’1  -L  “2  ^L  ^3  -L  ’  (A-1) 

where  u)  is  the  angular  velocity  of  a  launcher-fixed  coordinate  frame  and 

/V  /x  /\ 

1l»  and  form  a  set  of  dextral,  orthogonal,  unit  vectors  directed  along 
the  coordinate  axes.  The  matrix  counterpart  of  Eq.  (A-1)  is 

T 

w  =  (cOj^  W2  u)^)  ,  (A-2) 

where  the  superscript  T  denotes  the  matrix  transpose.  Of  course,  if  a 
matrix  is  to  be  used  in  place  of  Eq.  (A-1)  some  basis  must  be  adopted. 

In  the  case  of  Eq.  (A-2),  it  is  the  launcher-fixed  basis.  Herein,  the 
basis  adopted  for  writing  the  matrix  equivalent  of  a  vector,  or  dyadic, 

f)2 


f)3 

should  be  apparent,  in  some  cases,  from  the  context  of  the  equation  in 
which  it  appears.  In  other  cases,  the  basis  is  defined  in  the  text;  and, 
in  still  others,  the  basis  used  is  indicated  by  the  letter  L  or  R  written 
underneath  a  matrix  to  denote  launcher -fixed  (1,),  or  rocket-fixed  (R)  basis. 

Dyadics,  such  as  the  Inertia  dyadic  of  the  launcher  are  indicated  by 
double  underlines.  For  example,  the  inertia  dyadic  of  the  launcher 
referenced  to  the  launcher  support  point  0  is 


where  the  I^j  are  moments  (i=j)  and  the  negatives  of  products  (it^j)  of 
inertia  of  the  launcher  about  0. 

The  matrir;  equivalet^t  of  the  inertia  dyadic  1^^  is  the  Inertia  matrix. 


(A-4) 


where  it  is  obvious  that  the  launcher-fixed  basis  has  been  adopted. 


The  vector  cross  product,  wxR,  where  R 
matrix  equivalent  (d  R,  where 


^1  1,  ii,  the 

‘'C  ^ 


(A-5) 
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T 

and  R  =  (x  y  z  )  . 

If  It  is  convenient  to  use  certain  bases  to  write  the  matrix  equi¬ 
valents  of  vectors,  dyadics  and  vector  products,  then  the  matrix  forms 
must  be  transformed  before  matrix  products  can  be  formed.  For  example, 
let  A  denote  the  direction  cosine  matrix  which  defines  the  orientation  of 
the  rocket-fixed  system  Cxyz  with  unit  vectors  i,  j  and  k  and  let  u  denote 
the  matrix  equivalent  of  a  rocket-fixed  (with  respect  to  rotation)  vector, 

then  the  matrix  product  uj  A  u  and  the  vector  product  wxu  are  equivalent. 

T  ' 

Furthermore,  the  matrix  product  A  u  ^  w  is  equivalent  to  the  vector 
product  uxw.  The  basis  of  the  result  in  both  examples  is  the  launcher- 
fixed  one. 

In  transforming  vector/dyadic  equations  to  matrix  equations,  care 
must  be  taken  to  make  sure  that  two  quantities  with  different  bases  are 
not  multiplied,  added  or  subtracted.  To  illustrate  this,  we  consider  the 
velocity  of  a  point  P  in  the  rocket  (see  Fig.  Al) .  In  Fig.  Al,  0  is  a 
fixed  point  and  C  and  C  are  centers  of  mass  of  the  launcher  and  rocket, 

Lj 

respectively. 


Fig.  Al  Definition  of  vectors. 
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We  let  R  and  s  be  the  launcher-fixed  basis  equivalents  of  the  vectors  R  and 
—  —a  — 

L  L 

s  .  We  also  let  u  and  p  denote  the  rocket-fixed  matrix  equivalents  of  the 
R  R 

vectors  uandp.  We  also  use  to  denote  the  rocket-fixed  basis  matrix 
equivalent  of  the  rocket's  angular  velocity.  Then  in  matrix  form,  the 
position  of  P  (referenced  to  0)  is  given  by 


r„=R  +  s  +U  +  P 
-P  -  -a  -  - 

and  the  velocity  of  P  is 


(A-6) 


Vp  =  wx(R4«^)  +  (ds^/dt)^  (A-7) 

+  nx(u+p)  +  [(du/dt)j^  +  (dp/dt)p]  , 

where  (d*/dt)  and  (d*/dt)  denote  the  time-rates-of-change  of  the  launcher- 
L  K 

fixed  and  rocket-fixed  components  of  *;  i.e.,  "local"  derivatives. 

A  matrix  equivalent  of  (A-7)  is 

Vp  =  (ji)(R+s  )  +  a"^  Q(u+p)  +  g  +  A^  u  +  a"^  £  >  (A-8) 

where  x  and  x  denote  the  matrix  equivalents  of  the  local  derivatives 
L  R 
defined  above. 

Certain  special  matrices  are  used  in  this  appendix.  Some  are  used 
to  extract  needed  information  from  a  matrix  equation  in  order  to  remove 
constraint  forces  and  torques.  By  definition. 


1 


il 


0 


0 

0 


0 

0 


(A-9a) 


,0 


0 


0 


66 


and 


123 


0  0  0 
0  10 
0  0  1 


These  are  recognized  as  subsets  of  the  identity  matrix. 


E  = 


10  0 
0  1  0 
0  0  1 


(A-9b) 


(A-10) 


Clearly,  Ej^  can  be  used  to  retain  only  the  first  element  of  a  3x1  matrix, 
while  E^^  can  be  used  to  retain  the  last  two  elements  of  such  a  matrix. 


A. 3  Physical  Models 

Launcher/rocket  systems  are  relatively  complex  dynamic  systems.  A 
typical  system  consists  of  a  vehicle,  the  launcher  proper  and  one  or  more 
rockets.  The  vehicle  types  include  massive  track  vehicles,  trucks,  towed 
carriages  and  human  beings.  The  launchers  range  from  heavy  pods,  with 
mechanisms  to  orient  them,  to  simple  hand-held  tubes.  The  rockets  may  be 
of  several  types  as  far  as  size  and  geometric  configuration  are  concerned, 
but  are  generally  propelled  by  solid  propellant  rocket  motors.  A  physical 
model  which  incorporates  the  salient  characteristics  of  all  possible 
systems;  l.e.,  a  generic  model;  would  be  tremendously  complex.  For  this 
reason,  it  appears  best  to  adopt  physical  models  which  are  most  suitable 
for  the  investigation  of  the  particular  problem  of  Interest,  or  which  are 
dictated  by  a  specific  system  concept.  Here,  we  are  interested  in  de¬ 
termining  the  degree  to  which  a  launcher  can  be  expected  to  act  as  a 
passive  control  device  which  reduces  the  flight  errors  caused  by  imperfec¬ 
tions  of  the  rocket (s)  launclicd  from  it.  Relatively  simple  models  arc 


sufficient  for  tills  purpose. 


Launcher  Model 


The  physical  model  adopted  for  the  launcher  is  a  single  rigid  body 
with  two  degrees  of  rotational  freedom  (see  Fig.  A2).  Although  simple, 
this  launcher  model  is  adequate  for  light-weight  launchers  and  for  the 
sub-launcher  concept  described  in  the  body  of  this  report.  Only  single¬ 
round  launchers  are  considered. 

Rocket  Model 

The  rocket  model  is  a  variable-mass  body,  composed  of  a  rigid 
shell,  a  rigid  warhead,  solid  propellant  and  the  fluid  part  of  burned 
propellant  still  within  the  shell  (see  Fig.  A3).  The  aerodynamic  charac¬ 
teristics  are  "typical."  The  principal  effects  of  time-varying  mass  are 
that  (1)  the  center  of  mass  of  the  rocket  moves  with  respect  to  the  rigid 
body  part,  (2)  the  mass  and  moments  of  inertia  of  the  rocket  change  with 
time  ind  (3)  the  flow  of  fluid  within  the  rocket  produces  jet-damping. 
Interface  Models 

The  interface  between  launcher  and  rocket  is  rigid.  Hence,  the 
rocket  pitches  and  yaws  with  the  launcher  during  the  "detent"  and  "guidance" 
phases.  Two  methods  for  imparting  spin  to  the  rocket  are  incorporated. 
First,  spin  may  be  generated  by  a  rocket-fixed  torque  such  as  that 
delivered  by  spin  vanes.  Second,  and  of  most  application  in  this  report, 
spin  may  be  Imparted  by  helical  rails  which  cause  the  rocket  to  rotate 
through  a  specific  angle  during  the  guidance  phase  of  flight. 

Both  non-tip-off  launches,  in  which  the  rocket  becomes  totally  free 
of  the  launcher  at  one  Instant  of  time,  and  tip-off  launches,  in  which 
forward  support  of  the  rocket  is  removed  at  one  time  and  at  a  later  time 


Fig.  A2  Launcher  physical  model. 


Fig.  A3  Rocket  physical  model. 
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the  rear  support  is  removed,  are  modeled.  The  guidance  phase  ends  at  the 
beginning  of  the  tip-off  phase  (see  Fig.  A4). 

A. 4  Equations  of  Motion 

The  present  derivation  of  equatic  ns  of  motion  follows  closely  that  of 
Ref.  6.  The  principal  differences  are  (1)  the  present  derivation  is  more 
brief,  (2)  the  use  of  helical  rails  to  impart  spin  was  not  considered  then 
but  is  now,  and  (3)  some  errors  in  the  results  given  in  Ref.  6  have  been 
corrected  here.  As  indicated  above,  the  total  time  from  ignition  of  the 
rocket  until  its  Impact  is  divided  into  four  phases:  detent,  guidance,  tip- 
off  and  free-flight.  The  detent  phase  begins  prior  to  ignition  and  ends 
when  the  mechanism  which  constrains  the  rocket's  translation  along  the 
launcher  releases.  The  guidance  phase  begins  when  the  detent  phase  ends  and 
ends  when  Che  forward  support  of  the  rocket  releases;  l.e.,  tip-off  begins. 
The  tip-off  phase  ends  when  the  rocket  is  no  longer  in  mechanical  contact 
with  the  launcher. 

Because  the  number  of  degrees  of  freedom  varies  from  phase  to  phase, 
each  is  considered  separately.  However,  most  of  the  vector  notation  and 
the  coordinate  systems  utilized  are  shown  in  Fig.  A5.  The  OXYZ  system  is 
an  Inertial  one  and  has  its  origin  located  at  the  pivot  point  of  the 
launcher.  The  launcher  has  two  degrees  of  freedom  in  rotation,  yaw 
(angle  Ot  )  and  pitch  (angle  The  launcher— fixed  system  C  x  y  z  has 

J  ^  Li  1j  1j  Li 

its  origin  at  the  launcher  mass  center  and  its  x^^-axis  directed  parallel 
to  the  direction  of  travel  of  the  rocket.  The  z^^-axis  lies  in  a  vertical 
plane.  Vector  R  locates  Cj^  with  respect  to  0.  The  rocket-fixed  ax'y'z' 
system  has  its  origin  a  on  the  geometric  centerline  of  the  rocket  at  the  aft 
support  location.  Vector  s  locates  point  a  relative  to  C  . 

“3  L 


The  x'-axis 


r 
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Fi^.  A5  Notation  and  reference  frames. 
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lies  along  the  geometric  centerline  of  the  rocket.  The  rocket's  center  of 

mass  is  C  and  the  Cxyz  reference  frame  rotates  with  the  rocket.  Because 

the  mass  of  the  rocket  varies  with  time,  C  translates  with  respect  to  a. 

However,  the  Cxyz  system  remains  aligned  with  the  ax'y'z'  system  at  all 

times.  The  vector  u  locates  C  relative  to  point  a. 

The  vectors  >  ,  p  ,  (■  ,,  1’,,  and  j.;  are  shown  in  Fig.  A5.  A  point  on 
- 1  -m  -M  - 1  - 

the  line-of-act  ion  of  F,,,  is  located  by  8,  .  Tl.e  vectors  p  and  p  are  generic 
position  vectors  which  locate  elements  of  mass  dm  and  dM,  respectively,  in 
rocket  and  launcher.  Finally,  12  and  ui  are  the  angular  velocities  of  the 
rocket  and  launcher,  respectively. 

A. 4. 2  Detent  Phase 

During  the  detent  phase,  the  launcher  and  rocket  move  as  one  body, 
except  the  rocket  may  rotate  about  the  x'-axis.  If  we  let  Ij^^q  1 
denote  the  inertia  dyadics  of  the  launcher  about  0  and  the  rocket  about  C, 
respectively,  the  angular  momentum  of  the  system  about  0  is 

=  1^/0*-  i’-  ’  (A-11) 

where  r  =  R  +  s  +  u.  Because  (du/dt)  and  (d^u/dt^)  should  be  small  during 
“  •“  “3  ~  R  *“  R 

the  brief  detent  phase,  and  since  (ds  /dt)  F  (d‘s  /dt^)  H  0  during  this 

— <1  L  ”3  L 

phase,  the  equation  of  rotational  motion  of  the  system  about  0  is 

Tq  =  Il/o*  +  w  X  +  1  •  (dl2/dt)j^ +5  X  I*Sl  +  rxmr  ,  (A-12) 

where 

if  =  (dw/dt)  x(R+s  )  +  wx((iix(R+s  )l+(dS2/dt)  +  12  x(12xu)  .  (A-13) 

—  —  L  ^  —  —  — 

The  torque  T^  includes  a  constraint  torque  which  prevents  launcher  roll. 

Also,  in  matrix  form,  rocket  basis. 


Q  —  A  +  A  u)  , 


where  X  =  (p  0  0)  is  the  spin  rate  of  the  rocket  and 


0  0 


A  =  I  0  s0^  . 

LO  -sU^ 


(A-14) 


(A-15) 


Here,  we  have  used  the  shorthand  notati.^n,  r()^  =  rosO^^  and  s0^  E  siny^. 
Also, 


01  =  p  • 

From  (A-14)  and  (A-15)  the  matrix  counterpart  of  (df)/dt)  is 

R 

•  •  •  *' 

=  A+Aoj  -  AAo) 

R 

and  the  equation  of  rotational  motion  for  the  rocket  is 


I  .d^/dt)  +  5x|*9  =  T  .  , 


(A-16) 


(A-17) 


(A-18) 


where  T  consists  of  constraint  torque,!  ,  and  spin  torque,  T  ,  if  present 


The  constraint  torques  are  orthogonal  to  the  x-axis,  so  that 


T  =  E  {I  n  I  . 

R  "  " 


(A-19) 


Also,  since  A  =  (p  0  0)  , 


A  =  E,  {-A  (L  +  A  A  w  -  r'E^,lI  fi  +  Sl  1  fi]  +  l”  T  } 
-  =1  =  _  -  =  -  =  =23  =  -  -  =  -  =  -s 


(A-20) 


Letting  1  denote  the  moment  of  inertia  of  the  rocket  about  the  x'-a.cis 
s 

and  defining 


=  E  -  li  f'h 
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wc  have 

^  =  T  [X  A  u)  -  A  Jj)  -j"*E,r'E,S2Il7  +  (T  0  0)^}  .  (A-22) 

Ig  =1=  "  =  ~  =  ■  =  =1“  =l-=-  s 

The  translation  of  C  is  governed  by  the  equation, 

mr  =  ,  (A-23) 

GX  C 

where  F  is  the  external  force  on  the  rocket  and  F  is  a  constraint  force  which 
includes  the  detent  force.  In  matrix  form,  launcher  basis. 


F  =  -F  +m{-ru)  +  oxjCR  +  s  ) 

L  L  -  _  —  a 

+  ^  u  -  A^u  X  +  A^  u  X  A  u)}  ,  (A-24) 


~  T  ~ 

where  r=R+s^+A  uA. 

By  using  (A-22)  in  the  matrix  form  of  (A-12),  we  obtain  the  result. 


•  .  .j.  .. 

<u  =  -  01  <^1  +  A  (I  X  A  U)  -  Q  1  i)] 


-m  r  w  (o(R+s  )  -m  if  A  (uXM)  +  i^] 


(l/I  )[A  I  -mfAuKE  (IXAw  -  J'^E^T'e 


+  (T  0  0)}+  T,  , 

s  -0 


(A-25) 


where  the  "inertia  matrix"  is  given  by 


L,  =  I,  -mrr  +  A  I{E-(1/I  )I  E  )A 
=1  =L/0  —  =  =  =  s  =  =1  = 


(A-26) 


The  torque  T^  is  due  to  gravity,  rocket  thrust,  and  the  spin  torque  T^.  It 
also  Includes  a  constraint  torque  which  prevents  rolling  of  the  launcher. 
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From  the  geometry  of  Fig.  A6,  with  ^  0, 


and 


+  02  sina^  , 


u)2  =  02  cosa^  cosa^  +  sina^^ 


=  -a^coso^  slna^  +  o^  cosa^. 


If  =  0, 


(jjj  =  sina^  , 


=  ot^  cosa^ 


and 


(O3  =  03  . 


Hence , 

0J3  =  (JO2  tana3 

and 

•  •  /  2 

0)^  =  0J2  tana3  +  OJ2OJ3/COS  03 


(A-27) 


(A-28) 


(A-29) 


•  • 

By  substituting  (A-28)  and  (A-29)  into  (A-25)  one  may  solve  for  (^2  and  0)3. 
Then,  p  can  be  found  from  (A-22) . 

The  part  of  Tq  due  to  "external"  forces  is  T^  .  We  assume  that  the 
launcher's  rotation  is  restrained  by  linear  springs  and  viscous  dampers. 
The  torque  due  to  these  is 


It’D  °  ’  (A-30) 

where  K  and  C  are  constant  stiffness  and  damping  matrices,  respectively,  and 
T 

a  *  (aj^Aa2  0I3)  .  Here,  A«2  is  the  difference  in  n,  and  the  quadrant  elevation 


(QE)  angle. 
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The  torque  about  0  due  to  gravity  Is 


T  =  mrBg+MRBg, 

s 


(A-31) 


where 


1  0  0  1  f  sa3  “ir  ca^  0  -^^2 


B  =  lO  ^“3  ^11^  ^  ^ 

Lo  -sa^^  ca^J  L  0  0  iJ  L  sa^  0  ca^ 


and  g  =  (0  0  g)  .  Note  that  is  not  assumed  to  be  identically  zero,  but 

is  constrained  to  be  zero  for  this  investigation. 

Action  of  the  thrust  on  the  rocket  produces  a  torque  about  0  which  is 
given  by 


It  *  ?T  • 


(A-32) 


Finally, 


—  T  +  T  +T 

-0  -SD  -g  -T  -s  ' 


(A-33) 


The  kinematic  equations  for  the  launcher  for  this  and  subsequent  phases 
are,  if  E  0, 


(A-34a) 


«!  =  0)3  . 


(A- 34b) 


However,  the  computer  code  in  which  the  equations  are  incorporated  contains 
the  more  general  equations. 


^2  ^**^2  ”  ‘'*3  sin(ij^) /cosrt^  , 


(A-35a) 


=  (1)2  sinoij^  +  cosaj^ 


(A-35b) 
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and 


sina^ 


(A-35c) 


A. 4. 3  Guidance  Phase  Equations 

During  the  guidance  phase,  allowance  is  made  for  spin  due  to  helical 
rails,  or  the  spin  torque  T^.  Two  sets  of  equations  are  needed  because  the 
constraints  present  when  spin  is  by  helical  rail  are  different  from  those 
which  exist  when  there  is  a  rocket-fixed  spin  torque  due  to  spin  motor, 
spin  vanes,  or  a  similar  device. 


Spin  Torque  Case 

When  spin  is  generated  by  T^,  the  constraints  during  guidance  include 

a  force  =  (0  “  ^^1^  ^2^  ^”3^)^  denote  the  resultant 

external  force  on  the  rocket.  Then,  because  there  can  be  no  motion  in  the 

y  -  or  z  -direction, 

L  L 

=  m  E.,-{u)  w(R+s  )  -  (R+s  )(i) 

=ZJ  -  -  -  -a  --a- 

+  2  w  -  A^[u  n  +  n  uj}  ,  ^ 


,ex 


+  f: 


+  n 


T 

where  s  =  (x  y  z  )  and  only  x  varies  with  time.  Equation  (A-35)  can 

*"3  Lt  Lt  Li  1> 

a  a  a  a 

be  used  along  with  (A-22)  and  (A-23)  to  obtain  the  following  equation  for  w: 
L2  (L  =  -  w  u)  -(1/1^)  lA^  1  -mr  |23A^uH  E^I I AA  w 


p  I  fill  + 

"t‘ 

s 

0 

+  f 

0 

.0. 

.0 

T  ' 

+  A  [I  AAaj-12I12|  -mr  E-.l(iJw(R+s  ) 
=  =-=-  -=Z3  —  --a 


(A-36) 
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+  2  Hi  s  +A'^(uAAa)  +  fii7u)] 

+  (A-36  cont.) 

-0  ’ 


where 

L,  =  I,  +  a"^  I  A  -mi  E„„(f  -(l/I  )A^  G  E,  1  A] 

=2  =L/0  =  =  =  -  =23  -  s  =  -  =1  =  = 

-d/I  )a'^  I  E,  I  A  . 

Translation  of  the  rocket  is  governed  by  the  equation, 
T  ex  T 

m(x_  0  0)=  (F  0  0)  E,  {axo(R+s  ) 

L  1  “"i  ~  “*a 

a 

T*r* 

-(R+s  )(jL)  +  2  (jj  s  -  A  u[A  -  X  A  a) 

-  -a  -  -  -a  ss  -  -  -.  =  -. 

+  A  w]  +  a"^  §  ^  u  }  . 


(A- 37) 


(A- 38) 


Rotation  of  the  rocket  relative  to  the  launcher  may  be  found  from 
(A-22)  and  9^^  =  p. 

Helical  Rail  Case 

The  helical  rail  is  assumed  to  have  a  constant  pitch  such  that  if  the 
rocket  translates  a  distance  it  rotates  through  an  angle  l  C  about  its 
x'-axis.  It  follows  that  the  rocket  angular  velocity  relative  to  the 
launcher  is  (launcher  basis  or  rocket  basis) 

X  =  (c  x^  0  0)^  .  (A-39) 

a 

A  portion  of  the  helical  rail  is  shown  in  Fig.  A7.  The  angle  a 

H 

is  defined  by  the  following  equations: 

cosoj^  =  e(d/2)/[l+(Ld/2)n‘*  (A-40a) 

8ina„  =  l/[l+(cd/2)'’l^ 

H 


(A-40b) 
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Also  shown  in  Fig.  A7  are  the  F  ,  F  ,  F  and  F  ,  the  normal,  tangential. 

axial  and  spin  components,  respectively,  of  the  force  on  the  rocket  due 

to  the  rail.  Only  one  rail  is  shown;  however,  in  most  launchers  there 

will  normally  be  two  or  more.  Herein,  we  assume  two  diametrically  opposed 

rails,  so  that  the  total  force  on  the  rocket  due  to  the  rails  is 

2(F  +  F  )  =  2(F  +  F  )  . 

-n  -t  -a  -s 

From  the  geometry  of  Fig.  A7 ,  we  have 


F  =  -F  slna,.  -F  cosot,,  (A-41a) 

a  t  H  n  H 

and 

F  =  -F  cosa,,  +  F  sincx,,.  (A-41b) 

s  t  H  n  H 


We  assume  that  F  =  M,.F  ,  where  p.  is  the  coefficient  of  dynamic  friction 
t  f  n  f 

corresponding  to  the  rail/shoe  combination.  It  then  follows  that  the 
torque  on  the  rocket  due  to  the  rail  is  given  by 

^rail  =  '  ^^f  ’  <^-42) 

where  d  is  the  rocket's  diameter  at  the  shoe  location.  Equation  (A-42) 

may  be  used  to  eliminate  F  from  (A-41b)  and  obtain  the  axial  force  on  the 

n 

rocket , 


F  ,  =  -K  T  ,,  , 

axial  r  rail 

where 

=  (2/d)  (coscxj^  +  Pj  sin  jj)/(sinaj^  -  p^  cosot^^)  . 


(A-43a) 

(A-43b) 


The  equation  of  rotational  motion  of  the  system  about  0  is  again 

(A-12)  except  that  r  has  the  more  general  form  implied  by  (A-38);  l.e., 

s  is  not  constant, 
ra 
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Rotation  of  the  rocket  is  governed  by  the  equation, 

I  •  (df2/dt)„  +  fJxI  •  +  T  (A-44) 

=  “  R-=-  -  -  -rail 

where,  as  before,  T®*  is  the  external  torque  about  C  and  is  the  con¬ 
straint  torque  which  forces  the  rocket  to  pitch  and  yaw  with  the  launcher. 
Furthermore,  the  equation  which  determines  the  motion  of  the  rocket's 
center  of  mass  is 

mr  =  F®*  +  f‘^  +  F  .  ,  ,  (A-45) 

-  ~  -  -axial 

©X  c 

where  F  =  F_  +  F  and  F  is  a  force  which  constrains  the  motion  to  be 
-T  -g 

in  the  y.  -direction. 

L> 

By  writing  (A-44)  in  matrix  form  (rocket  basis)  and  solving  for 
we  find  that 

T  =  E  fI(A  w  -  X  A  0)  +  X)  +  51  1  fi]  .  (A-46) 

-tail  =1  -  =  L - -  -  =  - 

This  result  for  T  can  be  used  in  the  matrix  counterpart  of  (A-45)  to 
-rail 

get  the  result, 

••  — 1  -w*  T~  * 

s  =  A  E,  {-2u)s  -  (jjciis  +  A  u  A  o) 

-a  =1  — a  - a  =  "  =  r 

li 

-A^(uX^  +  5151u  -  F.j,/m)  -  Bg 

-Ej.  [I (Aw  -  XAw)  +  Ql^]}  ,  (A-47) 

where  A=1  -e(l  0  0)A^  u  A(1  0  0)^  +  cK  1  /m. 

as  •  ss  i;-  S 

The  expression  (A-47)  can  be  used  in  the  matrix  counterpart  of  (A-12) 


to  get 
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•  •J* 

L  a)=-djl,  --(jJ-A  +  A  I  A  A  cj 

=rail  —  -  =L/0  -  =  -  =  =  -  =  - 

-[mr(E  -  -c(A^  1  A  -mr  A^  u  A)A"^E^]{2w 

'J*  -N,  y  ~~ 

+  ioajr+()iA  Au+A  AAu} 

-[mr  +  c(A^  I  A  -mr  a"*^  ]A~‘  ^ 

+  Tq  .  (A-48) 


Here, 

Wall  *  Il/o  +  4’'  1  i  -  "££ 

+  A  *  [mr  +  g(a”^  1  a  -  mf  a"^  u  A)  ]E^  r 

-  A"^  r  I  A  .  (A-49) 


Then,  p  =  e  x  . 

a 

A. 4. 4  Tip-Off  Phase  Equations 

During  the  tip-off  phase,  if  any,  the  rocket  is  assumed  to  rotate 

with  respect  to  the  launcher  about  the  point  a.  The  procedure  for  finding 

the  equations  of  motion  for  this  phase  is  the  same  as  for  the  previous 

phases  except  that,  since  there  is  no  constraint  torque  about  point  a 

(friction  neglected)  the  derivative  (d  Q/dt)  can  be  computed  directly. 

The  main  result,  from  which  <v  can  be  found  is 

L 


L,  0) 

=3  L 


-  u)  w  -(R  +  s^)E*  m{2  w  s^ 

+  axii(R+s^)  +  A^ilfiu  -  Bg  -  A 


(A- 50) 


(A-50,  Cont.) 


r 


H4 

+  (R+s  )E„.  m  u  1  +(?^-u)F^  +  T  } 

■"  ““Si  ^  “•  “  ■“cl  “X  “  ”X  ”S 

+  lo  -  “■  18  -  . 


where 


^  -'"(K+.s  )i:*(K+s  ) 

=3  =L/0  -  -a  =  -  -a 

E*  =  E  +  E  A^  mu  I  *  u  A  E  , 
=23  =23  =  -  =a  -  =  =23 

T 

I  =  I  -mu  A  E  A  u  , 

=a  =  -  =  =23  = 

T 

A  =  C  B  , 


"1 

0 

0" 

'cO 

0 

-s()“ 

■  cT 

sT 

O' 

c  = 

0 

C'l' 

S<t) 

0 

1 

0 

-sT 
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.0 

-s'! 
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.  sO 

0 

cO. 

_  0 

0 

1. 

(A-51a) 

(A-51b) 

(A-51c) 

(A-Sld) 

(A-51e) 


and  M  is  the  launcher  mass.  Tlie  Euler  angles  appearing  in  C  are  the  yaw 
(4'),  pitch  (0)  and  roll  (O)  angles  of  the  rocket  with  respect  to  the 

fixed  reference  frame  OXYZ  (see  I'ig.  A8) . 

• 

The  equation  for  U  is 
R 

=  l”'{-  [2  il  +(SL^,  -  u  A  E^ 

R 

+  T  +  m  u  A  E  [Bg  +  (R+s  )a) 

-s  -  =  =23  =-  -  -a  v 

-  U)a)(R+s  )  -2  to  s  )  }  ; 

—  -  -a  -  -a 


(A-52) 


